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We present self- consistent calculations of electron transport in graphene nanoconstrictions within 
the Hartree approximation. We consider suspended armchair ribbons with V-shaped constrictions 
having perfect armchair or zigzag edges as well as mesoscopically smooth but atomically stepped con- 
strictions with cosine profiles. Our calculations are based on a tight-binding model of the graphene 
and account for electron-electron interactions in both the constriction and the semi-infinite leads 
explicitly. We find that electron interactions result in (i) Electrons accumulating along edges of 
the uniform ribbon and along zigzag and cosine constriction edges but not along armchair constric- 
tion edges, (ii) The first subband showing almost perfect transmittance due to localization at the 
uniform graphene boundary except at low energies for the cases of zigzag and cosine constrictions 
where Bloch stop-bands form in related periodic structures, (iii) The second subband being al- 
most perfectly blocked by the constriction, (iv) Electron interactions favor intra-subband scattering 
while the non-interacting electron theory predicts predominance of inter-subband scattering, (v) 
Conductance quantization for the first few conductance steps being more pronounced for armchair 
constrictions but less so for zigzag constrictions, (vi) A much more prominent 2e 2 /h conductance 
plateau for the cosine constriction than is found in the absence of electron interactions. Possible 
implications for recent experiments are briefly discussed. 



I. INTRODUCTION 

Graphene nanoribbons are strips of graphene several 
nanometers wide and of arbitrary length. Their unique 
electronic structure and transport properties that arise 
from the linear, massless Dirac-like spectrum of the un- 
derlying honeycomb lattice of graphene are attracting a 
great deal of interest at the present timepJ In experimen- 
tal studies of graphene nanoribbons the electronic charge 
density in the ribbon is usually varied by the application 
of a variable voltage to a gate electrode located near the 
ribbonPHni if ? due to the application of the gate volt- 
age, there is a net charge on the ribbon, the electronic 
charge density in the ribbon ceases to be uniform, and 
there is a strong redistribution of the charge towards the 
edges of the ribbon. This charge redistribution and its ef- 
fects on the electronic structure and transport have been 
examined theoretically by several authors for ribbons of 
uniform width separated from the gate electrode by a di- 
electric film! 11 " 14 ! It has been predicted that the charge 
redistribution results in a 1/ y/x charge singularity^ at 
the edge of the ribbon (here x is the coordinate normal 
to the ribbon bounda ry) and modification of the elec- 
tron dispersion relation.^ Electron transport in graphene 
nanoribbons with constrictions (GNCs) is also attract- 
ing theoretical^ and experimental attention at the 
present time. However, the effects of charge redistribu- 
tion in GNCs have not as yet been discussed in the liter- 
ature. In this paper we explore this topic theoretically. 

The present theory accounts for the effects of the 
Coulomb repulsion that gives rise to the charge redistri- 
bution within the self-consistent Hartree approximation 
that has been used previously to study charge redistri- 
bution in uniform ribbons .EHUl We consider infinite rib- 
bons with armchair edges and V-shaped or cosine shaped 
constrictions and treat the effects of electron-electron in- 



teractions at the Hartree level throughout these entire 
structures. The ribbons that we consider are suspended 
above a dielectric layer which covers a gate electrode as 
in a recent experiment.^ We treat the effects of the di- 
electric and gate within an image charge model. We 
describe the charged nanoribbon by fixing the chemical 
potential, and obtaining a self-consistent solution where 
all electronic states with lower energies than the Fermi 
energy are filled. From this calculation we obtain the 
Hartree electronic energy bands and charge densities, and 
compute the conductance. Our model considers only p z - 
orbitals of the graphene within the tight-binding approx- 
imation. The structures studied below include more than 
10000 carbon atoms inside the computational area. 

We find that strong charge accumulation along the con- 
striction boundary occurs or does not occur depending on 
the type of constriction that is considered. As more elec- 
tron subbands become populated the lowest subbands 
gradually localize near the graphene boundary. We find 
that, in the Hartree model, electrons in the first subband 
of the ribbon are transmitted almost perfectly through 
the constriction with little intersubband scattering, ex- 
cept at the small Fermi energies where they can be res- 
onantly reflected by a constiction having zizgag edges. 
By contrast, in the non-interacting electron model strong 
intersubband scattering occurs at the constriction. De- 
pending on the type of constriction conductance quan- 
tization is predicted to be more or less pronounced in 
the Hartree approximation than in the non-interacting 
approximation, a finding that may be relevant to recent 
experimental observations. 

The remainder of this paper is organized as follows: 
Our model and method of solution are described in Sec- 
tion [Hj Our numerical results are presented in Section 
Mil We summarize our main conclusions in Section HVl 
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II. MODEL 

We consider a graphene nanoribbon (GN) suspended 
in a way similar to that in the experimental setup of 
Tombros etalP^ The GN is separated from the back gate 
by layers of a dielectric and air, see Fig. 1(a). For the 
dielectric material we choose Si02 with relative permit- 
tivity e = 3.9. The NR is attached at its two ends to 
semiinfmite leads represented by ideal ribbons having the 
same width W as the NR. Four different types of devices 
are considered in the following: an ideal uniform rib- 
bon, and ribbons having V-shaped armchair and zigzag 
and constrictions or constrictions with cosine profiles im- 
posed on one side. We disregard any defects other than 
the atomic steps at the boundary of the cosine-shaped 
constriction in the present study. The host configuration 
is taken as armchair as are the edge configurations of the 
semiinfmite leads. 

As representative devices we consider GNs of width 
W = 10 nm and length L = 27 nm, Fig. 1(a). There 
are 92 carbon atoms in the cross-section making the rib- 
bon semiconducting. The constrictions, if imposed, are 
V-shaped or cosine-shaped trenches 5 nm deep inside the 
ribbon thus leaving half of the width W of the ribbon in 
the narrowest part of the constriction for electron propa- 
gation. The length L is taken long enough to include part 
of the leads near the constriction and therefore to treat 
the leads accurately. Because for the V-shaped constric- 
tions we consider atomically ideal boundaries, the apex 
angles of these constrictions are 60 and 120 degrees for 
armchair and zigzag cases, respectively. The Si02 and air 
layers are both 50 nm thick. Thus the back gate is 100 
nm from the nanoribbon. We performed simulations for 
different ribbon widths for these constriction geometries 
and all of the results showed similar features. 

The system shown in Fig. 1 is described by the Hamil- 
tonian 



(a\dj + h.c^j , (1) 
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FIG. 1: (a) Representative device with freestanding graphene 
nanoribbon. (b) Image charge model. 
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TABLE I: Coordinates and electron and image charge densi- 
ties for the model shown in Fig. [IJb). k — refers to the 
electron density no that gives rise to the image charge densi- 
ties for which k > 0. d is the distance between the graphene 
and the dielectric and / is the thickness of the dielectric. 



nearest-neighbor atoms; V/ 1 is the Hartree potential at 
atom i which results from the Coulomb interaction with 
the uncompensated charge density —en in the system 
(including the image charges). In coordinate space the 
Hartree potential can be written as 



V H (r) 



r-/|2 . 



A' 



(2) 



where — en/ c (r / ) is the k th electron or image charge placed 
at distance b^ from the graphene layer. The image 
charges included in the m odel keep the back-gate elec- 
trode at zero potential.^ ! 13 1 20 1 The potential due to a 
charge density —eno located a distance d above a dielec- 
tric (with dielectric constant e and thickness I) that is 
over a metal gate, as shown in Fig. [ljb), can be described 
by an infinite number of image charge densities. The 
first few image charge densities and their ^-coordinates 
(measured from the position of the electron charge) and 
also the electron charge density itself are given in Table 
[TJ The first row (k = 0) describes the direct Coulomb 
interaction between electrons in the graphene layer. Be- 
cause the contributions from the image charges decrease 
rapidly as k grows and in order to facilitate computation 
the results presented below were obtained keeping only 
the k = 0, 1, 2, 3 terms. 

The integration in Q was performed over the whole 
device including the semiinfmite leads. In order to in- 
clude electron-electron interactions over the whole sys- 
tem, we partition the system into three parts, the inter- 
nal computational region and two semi-infinite leads 1211221 
The internal region incorporates not only the constriction 
but also segments of uniform ribbon on both sides of it, 
including part of the leads. The semi-infinite leads them- 
selves begin far enough from the constriction to ensure 
that the total self-consistent potential and the electron 
density do not change appreciably along the leads, i.e., 
the electron density and the potential in the leads are 
not affected by the internal region. Thus the leads can 
be considered as uniform graphene ribbons. 

Starting from the Hamiltonian Eq. ([I]), we evaluate 
the Green's function numerically using the technique de- 
scribed by Xu et alPH The Green's function in the real- 
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space representation, G(r, r), provides information about 
the local density of states at site r, 

LDOS(r,E) = -^[5(r,r,E)] , (3) 

where factor 2 takes account of the spin degeneracy and 
S is the area corresponding to one carbon atom. The 
density of states is DOS(E) = J dr LDOS(r,E). The 
LDOS can be used to calculate the electronic density n(r) 
at site r 

poo 

n(r)= / dELDOS(r,E)f(E-E F ), (4) 
Jvc 

where Ep is Fermi energy and / is the Fermi-Dirac dis- 
tribution function. All the calculations reported in this 
paper correspond to the temperature T = 10 K. A zero 
temperature version of Eq. [4] has been used previously in 
Ref. [13j In general at non-zero temperatures the lower 
limit of integration in Eq. [4] is — oo. However, in the 
present work the LDOS is zero in a range of energies of 
width much larger than kT below the charge neutrality 
point V c and therefore Eq. [4] is a good approximation. 
The position of the charge neutrality point V c at a given 
Fermi energy is determined numerically from solution of 
the Schrodinger equation 

H^ = E(k)^, (5) 

with H being the Hamiltonian ([!]) and the wave function 
\I/ obeying the Bloch theorem 

* m +M = e ifcM * m , (6) 

where k is the Bloch wave vector and \I/ m is the Bloch 
wave function at coordinate m; M = 3a is unit cell length 
of the armchair ribbonP^HH] Having calculated the Bloch 
states and constructed the band diagram one can readily 
obtain the number of the Bloch states N Bloch for a given 
Fermi energy that in turn serves as a basis for analysis of 
transport properties of GNCspS Because the projections 
of the two Dirac points in the armchair ribbon coincide 
at k = we solve ([5| at zero wave vector and find V c 
from eigenvalues E(k = 0). 

The integration path in Q goes along the real axis and 
a fine integration grid is used to capture the locations 
of the subband edges and quasibound states if any are 
present. 

Since the Hartree potential Vh given by Eq. ^ de- 
pends on the electron density n(r) which is a solution 
of the Schrodinger equation with the Hamiltonian ([!]), 
these equations need to be solved iteratively. The iter- 
ation process is executed until the convergence criterion 
!!m < 10 -3 is met, where n 7 ^ and n™ ut are the in- 

put and output average values of the electron density at 
the m-th iteration. In the cases where a constriction is 
present in the GN the above computation proceeds in two 
stages. At the first stage, self-consistency is achieved for 
the uniform ribbon. Then constriction is imposed on the 
ribbon and self-consistency is achieved again. 



Having calculated the electron density and the position 
of the Dirac point numerically, we are in a position to find 
the conductance 

ij 

as a function of the Fermi energy. Here Tij(E) is the 
transmission coefficient from subband j in the left lead 
to the subband i in the right lead, at energy E. Tij(E) 
is calculated by the recursive Green's function method, 
see Ref. [23] for details. 

The Fermi energy and charge neutrality point are re- 
lated to a value of the gate voltage measured in experi- 
mental setup as V c + Ep — eV^! 11 * 13 ^ In a real device, it 
is the gate voltage V g that results in a change of carrier 
density in the graphene ribbon. We define eV g as the 
chemical potential difference between the metallic gate 
and the ribbon, necessary to accommodate extra carriers 
in the graphene and remove them from the metallic gate. 
Knowledge of both V c and Ep thus allows one to esti- 
mate the value of V g used in an experiment. Note that 
V c = in the non-interacting approach. 



III. RESULTS 

Figure 2 shows results for uniform GNs: (a) the con- 
ductance, (b) the charge neutrality point, (c) the po- 
sitions of the bottoms of the subbands that are near 
the Fermi energy, and (d) and (e) the wavevectors at 
the Fermi energy calculated within the Hartree and non- 
interacting models, respectively. The conductance values 
for the two models for the same values of the Fermi energy 
are very close to each other. In both cases the conduc- 
tance increases by the quantum 2e 2 /h each time a new 
subband opens for propagation. The opening of a new 
electron subband also results in an increase in the slope of 
the charge neutrality point V c vs. Fermi energy as is seen 
in Fig. |2|b); V c increases monotonically with increasing 
Fermi energy because the electron density on the ribbon 
increases. The slope change is caused by the additional 
contribution to the electrostatic potential on the ribbon 
due to the charges populating a new subband. Assum- 
ing a parabolic dispersion near the subband edge yields 
an E~ x l 2 divergence of the DOS that leads to an addi- 
tional electron density Sn oc J E~ x l 2 dE oc E 1 / 2 where 
E is measured from the subband edge. The associated 
charge contributes to the Hartree potential that in turn 
leads to the rise in V c . 

The subband energy position shown in Fig. 2(c) re- 
veals a linear drop of the energy levels relative to Ep 
as Ep increases. Each time an energy level crosses Ep 
electrons start populating the GN and contribute to the 
electrical conductance. Note that the bottoms of the 
subbands in the GN do not show any pinning to the 
Fermi level such as that observed in conventional quan- 
tum wires^ and open quantum dots^. As can be seen 
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in Fig. |2[d), the electron interactions modify the band 
structure of the GN: The results for the Hartree model 
show avoided crossings and the two lowest subbands hav- 
ing smaller velocities. The reason for these modifications 
of the band structure due to electron-electron interac- 
tions can be understood from the analysis of the charge 
distribution in the GN that is shown in Fig. [3j The elec- 
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FIG. 2: (color online) Results for uniform GNs: The conduc- 
tance (a), charge neutrality point (b), positions of the bot- 
toms of the subbands that are near the Fermi energy (c), and 
the wavevectors at the Fermi energy calculated within the 
non-interacting (d) and Hartree (e) approaches vs. the Fermi 
energy The inset in (c) shows the DOS for Ef — 0.105£; the 
peak marks the position of the bottom of the subband. The 
outset in (c) shows the derivative of the Fermi-Dirac distribu- 
tion function for T = 10 K. Arrows in (a) mark the energies 
used for charge density plots in Fig. [3] t — 2.7 eV. 
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FIG. 3: Comparison of the electron concentrations in uni- 
form GNs calculated in the non-interacting and Hartree ap- 
proaches for different numbers of populated subbands. The 
electron concentration in the armchair unit cell oscillates be- 
tween neighboring carbon atoms due to the specific structure 
of the wave functions, see e.g. Ref. [26j A half of the cross- 
section is shown. The electron Fermi energies for the different 
plots are indicated by arrows in Fig. [TJa). 



tron interactions in the Hartree model result in strong 
redistribution of the charges towards the edges of the rib- 
bon when the Fermi energy is increased. The larger the 
Fermi energy, the stronger redistribution of the electron 
density. Note that the charge accumulation along the 
boundaries of uniform graphene strips was also discussed 
in Refs. EHED 

Figure |4jc) shows the conductance as a function of the 
Fermi energy for a GN with an armchair constriction. 
The non-interacting electron calculation predicts conduc- 
tance quantization in steps of 2e 2 /h over the whole range 
of Ef values shown in Fig. Qc). However, the results of 
the Hartree calculation show only first two conductance 
plateaus to be well defined. The transition between these 
plateaus appears sharper than that predicted by the non- 
interacting calculation. In order to understand the reason 
for the better quantization for the first few conductance 
steps in the Hartree approach let's inspect the electron 
concentration distributions and transmission coefficients 
for Ef — 0.13t shown in Fig. [4] (a) and (b). The charge 
density is enhanced along the bottom straight boundary 
m Fig. [4| (a) and (b) in a similar way to that in the 
uniform GN discussed above. However the enhancement 
that occurs along the armchair constriction is very non- 
uniform: Much less charge accumulates near the apex at 
the narrowest part of the constriction; see Figs. Qa) and 
(b), respectively, for the electron concentration distribu- 
tion n H (r) in the Hatree calculation and the difference 
between the Hartree and non-interacting electron con- 
centrations n H (y) — n nom (r). The partial transmission 
Tij in the Hartree model reveals nearly perfect transmit- 
tance for the first state, see inset in Fig. |4^a). Note that 
Xii ~ 1 for all of the conductance steps in the Hartree 
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FIG. 4: (color online) (c) The conductance as a function of 
the Fermi energy for GN with armchair constriction: solid line 
- non-interacting electrons; line with dots - Hartree model, 
dashed line - (non-interacting) conductance without constric- 
tion that equals the number propagating states in the leads. 
The shaded gray area denotes the number, N Bloch of propa- 
gating Bloch states at the Fermi energy for the infinite mod- 
ulated ribbon in which the constriction is repeated period- 
ically, (a) The electron concentration n H calculated in the 
Hartree model, (b) The difference between Hartree and non- 
interacting electron concentrations n H - n noni . The Fermi 
energy in (a) and (b) is Ef — 0.13£ and is marked by the 
arrow in (c). The insets in (a) and (b) show the partial trans- 
mission probabilities Tij for the Hartree and non-interacting 
calculations, respectively. 



theory. However, the transmittance due to the second 
state was found to be strongly suppressed. for the 
non-interacting approach shows a fairly uniform distri- 
bution over all states, see inset in Fig. |4jb). To under- 
stand this phenomenon, let's consider the square moduli 
of the wave functions shown in Fig. [5] The first state 
in the Hartree approach transmits nearly perfectly be- 
cause of its localization near the straight bottom bound- 



ary. By contrast, the second state is mostly localized 
near the top boundary, where the constriction is located 
and transmission is therefore blocked. Note that local- 
ization develops gradually as Ep increases and the first 
and second states propagating along opposite boundaries 
reveal very similar dispersion for the uniform ribbon for 
Ep > 0.15t, see Fig. [2^d). These two states become 
mostly trapped within triangular wells at the straight 
ribbon boundaries. It is also worth noting that this phe- 
nomenon holds true for all of the constriction shapes 
studied below. It does not occur in the non-interacting 
electron model because of the absence of change accumu- 
lation near the boundaries. In general, we find electron 
interactions to favor intra-subband scattering whereas in 
the non-interacting approximation inter-subband scatter- 
ing predominates. The conductance step degradation in 
the Hartree model for Ep > 0.2t in Fig. [4jc) is related 
to the overall poor transmittance of the highest states, 
where the charge accumulation along the ribbon's edges 
but not in the constriction itself is further increased. 
When this happens the constriction becomes a more ef- 
fective obstacle to electron propagation through the rib- 
bon in the Hartree theory than in the non-interacting 
electron approximation. 
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FIG. 5: (color online) The wave function square modulus \^i\ 2 
for a GN with an armchair constriction calculated within the 
non-interacting and Hartree approximations, left and right 
columns, respectively. The Fermi energy is Ef = 0.13£, see 
the arrow in Fig. Etc). The corresponding partial transmis- 
sion probabilities for the zth-state are shown in the insets in 
Fig. |ja),(b). 

The conductance as a function of E F and a repre- 
sentative electron concentration distribution n H for a 
zigzag constriction are shown in Figures |6jb) and (a). 
The charge density along the zigzag constriction edge is 
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FIG. 6: (color online) (a) The electron concentration n H 
in the Hartree approximation. Concentrations equal to or 
greater than 8xl0 17 m -2 are colored yellow. The dashed 
rectangle bounds region with the electron concentration plot- 
ted on a different scale. The inset in (a) shows the partial 
transmission probabilities Tij in the Hartree model, (b) The 
conductance as a function of the Fermi energy for a GN with 
a zigzag constriction. The labels and meaning of the shading 
are the same as in Fig. [4|c). 



strongly enhanced by values up to an order of magnitude 
larger than for the armchair constriction; see the density 
scales in the inset in Fig. [6^a). One reason for this is the 
electron localization at zigzag edges previously predicted 
in both non- interacting and interacting electron theories; 
see, for example, Refs. I24H261 This is a topological prop- 
erty of zigzag-terminated ribbons. Another reason is the 
effect of electron interactions that increases the charge 
density along the edge further. The electrons occupy 
only one graphene sublattice along the zigzag edge while 
they occupy both sublattices along the armchair edge. 
We find the crossover between these charge occupations 
to occur over a distance of about 10 carbon atoms at 
the armchair-to-zigzag junction. The partial transmis- 
sion for the first subband is T\\ ~ 1 similarly to the case 
of the armchair constriction although it is suppressed at 
Ep ~ 0.05£ and ~ 0.08t due to resonant backscattering 
by strongly localized states at the zigzag edge. It is worth 
noting that in modulated ribbons consisting of periodi- 
cally repeated identical constrictions stop-bands form at 
these energies. These stop-bands, where the number of 
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FIG. 7: (color online) The same as Fig. [6] but for a GN with 
a constriction having a cosine profile. 



propagating Bloch states N Bloch (the grey shading in 
Fig. [6^b)) in the periodic structure is zero, are similar to 
those predicted in edge-corrugated graphene ribbons P^l 
However, here the stop-bands form due to electron in- 
teractions. The mismatch between the electronic struc- 
tures of the armchair host and zigzag constriction also 
contributes to less pronounced conductance quantization 
being seen in Fig. |6|b) for the zigzag constriction than 
in Fig. [4jc) for the armchair constriction. 

Figure [7] shows a representative electron concentration 
distribution n H and the conductance for a constriction 
with a cosine profile. We found strong accumulation 
of the charge density along the boundary of the con- 
striction itself, similar to that for the zigzag constric- 
tion. Note that even though the overall shape of the 
cosine constriction is smooth there are multiple pieces of 
the zigzag terminated edges connected to each other by 
atomic steps along the constriction boundary. The con- 
ductance dip at Ep ~ 0.08£ is correlated to the absence 
of propagating Bloch states at this energy in the cor- 
responding modulated ribbon with periodically repeated 
constrictions, as is the case for the similar features in Fig. 
6(b) for the zigzag constriction. Interestingly, although 
in Fig. [7|^b) only hints of conductance plateaus can be 
seen in the results for the non-interacting electron model 
(red solid curve), a very pronounced first conductance 
plateau is found when electron interactions are included 
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in the Hartree model (dots), although only hints of higher 
plateaus can be discerned in this case as well. We find 
the first conductance plateau to remain robust for GNCs 
with the cosine shape of length 7-16 nm (not shown). It 
may be relevant that in their experimental study of a 
constriction with an overall smooth profile (albeit for a 
much larger structure than those studied here) Tombros 
et al\^ observed a very pronounced first integer conduc- 
tance plateau (g ~ 2e 2 /h) but only very weak higher 
plateaus or only hints of higher plateaus. In such a sce- 
nario, a well defined first plateau may be caused by the 
lowest state being adiabatically transmitted along one of 
the graphene boundaries in an asymmetric constriction 
with curvature of one boundary being much smaller than 
that of the other. Because no image of the constriction on 
which the transport measurements were carried out was 
present ecP^, the degree of asymmetry is not known. Note 
also that the experimental device in Ref. [19] was fab- 
ricated by current annealing that substantially reduces 
edge disorder at both boundaries. 

IV. SUMMARY 

In conclusion, we have presented a self-consistent 
model of electron quantum transport in graphene ribbons 
and constrictions. The latter are represented as trenches 
of depth 5 nm and length range 7-17 nm and having dif- 
ferent shapes. The model is based on the Green's func- 
tion formalism and accounts for electron-electron inter- 
action within the Hartree approach. The Hartree model 
predicts several features not found in the non-interacting 
model. The electron charge density gradually accumu- 



lates along straight boundaries of uniform ribbons as the 
Fermi energy increases. However, accumulation at the 
constriction depends strongly on the details of the con- 
striction geometry. There is little if any charge accumu- 
lation along a V-shaped armchair constriction boundary 
but strong accumulation along a V-shaped zigzag bound- 
ary or a boundary with an overall smooth cosine profile. 
For each of these constriction types imposed on a ribbon 
with armchair boundaries, except near isolated reflection 
resonances, we find almost perfect transmission of elec- 
trons in the first subband (with little inter-subband scat- 
tering) and almost perfect reflection of electrons in the 
second subband of the ribbon in the Hartree model. By 
contrast, we find the constriction to induce strong inter- 
subband scattering of electrons for every subband in the 
non-interacting electron model. For the constriction with 
the cosine profile, the first integer conductance plateau is 
much more pronounced in the Hartree model than in the 
non-interacting model, a finding that may be relevant to 
the recent experiment of Tombros et a/pH The transport 
properties of two lowest subbands are the result of by 
electron localization near the opposite boundaries of the 
ribbon. An analogy can be drawn between perfect trans- 
mission along a uniform graphene boundary and the edge 
state transport in the quantum Hall effect, where both 
are immune to defects in interior of a device. 
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